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INTRODUCTION 


To  increase  the  maximum  pressure  a  cylinder  can  contain,  it  is  common 
practice  to  produce  a  more  advantageous  stress  distribution  involving  residual 
compressive  hoop  stresses  near  the  bore  by  autofrettage  treatment  of  the 
cylinder  prior  to  use  (ref  l).  The  determination  of  residual  stresses  is 
important  in  the  fracture  analysis  and  the  fatigue  life  estimation  (refs  2-6). 
There  is,  however,  considerable  disagreement  over  solutions  obtained  by 
different  investigators  for  the  residual  stress  distribution  in  the  cylinder 
after  the  autofrettage  process  (refs  7-12).  This  discrepancy  in  residual 
stress  is  a  result  of  different  mathematical  methods,  end  conditions,  and 
material  models.  Different  assumptions  for  the  material  properties  such  as 
compressibility,  yield  criterion,  flow  rule,  hardening  rule,  Bauschinger 
effect,  etc.,  can  lead  to  many  material  models.  Most  of  the  earlier  solutions 
for  residual  stresses  were  based  on  the  assumption  of  elastic  unloading  and 
only  a  few  considered  reverse  yielding  (refs  8,11).  For  unloading  with 
reverse  yielding,  there  is  no  general  agreement  in  the  literature  over  which 
material  model  should  be  used.  Many  plasticity  theories  have  been  proposed 
and  reviewed  (ref  13),  yet  no  theory  is  completely  adequate.  In  particular, 
it  seems  that  no  theoretical  model  has  been  given  to  represent  accurately  the 
actual  material  behavior  in  a  high  strength  steel  as  reported  by  Milligan, 

Koo,  and  Davidson  (ref  14). 


References  are  listed  at  the  end  of  this  report. 
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In  this  report  a  new  theoretical  model  Is  proposed  with  one  attempt  to 
give  a  better  representation  of  the  actual  loading/unloading  behavior  in  a 
high  strength  steel.  A  closed-form  solution  for  calculating  residual  stresses 
in  autof rettaged  thick-walled  cylinders  Is  obtained,  and  some  numerical 
results  will  show  the  Influence  of  the  Bauschtnger  and  hardening  effects. 

MATERIAL  BEHAVIOR  AND  MODELING 

The  material  chosen  for  this  Investigation  was  a  modified  4330  steel 
having  a  martensitic  structure.  A  description  of  Its  chemical  composition  and 
various  heat  treatments  is  given  In  Reference  14  by  Milligan,  Koo,  and 
Davidson.  They  studied  material  behavior  by  utilizing  a  uniaxial  tension- 
compression  specimen.  Figure  l  shows  the  stress-strain  curve  during  loading 
and  unloading  after  overstrains  in  tension.  The  stress-strain  curve  during 
loading  was  assumed  to  be  elastlc-perfectly  plastic.  This  was  a  good 
approximation  since  the  tensile  test  exhibited  very  little  strain-hardening 
during  loading.  This  would  also  be  true  for  other  steels  used  such  as  in 
References  6  and  9. 

Initially  the  yield  stresses  in  tension  and  compression  are  approximately 

equal  so  that  the  material  can  be  considered  as  isotropic.  However,  the  ratio 

of  the  yield  stress  upon  reverse  yielding  to  the  initial  yield  stress  is 

strongly  affected  by  overstrain  as  shown  in  Figure  l.  The  values  of  the 

Bauschinger  effect  factor  (BEF)  also  depend  on  the  offset,  and  the  0.1  percent 

offset  was  chosen  for  the  present  study.  Figure  2  shows  the  Bauschinger 

p 

effect  factor  (f)  as  a  function  of  percent  tensile  overstrain  (e  ).  The  graph 
shows  a  decrease  of  the  BEF  with  an  increasing  amount  of  tensile  prestrain  up 


to  approximately  two  percent,  at  which  point  it  becomes  effectively  constant 
(ref  14). 

The  Bauschlnger  effect  factor  is  very  important  in  determining  the  range 
of  elastic  unloading.  After  reverse  yielding  occurs,  a  very  large  slope  of 
strain-hardening  will  develop,  even  though  the  initial  tensile  test  exhibits 
very  little  strain-hardening.  A  bilinear  model  for  elastic-plastic  unloading 
is  proposed  here,  as  shown  in  Figure  1.  Choosing  a  new  coordinate  system 
(o',c')  with  the  origin  at  the  point  before  unloading,  we  have  for  the  plastic 
portion  of  the  reverse  yielding  curve 

o'/ o0  -  1  +  f  +  m'  C7(l-m')  (1) 

where  C'  -  (E/o0)e'P,  E  is  Voung' s  modulus,  0o  is  the  initial  yield  stress, 
tn'E  is  the  slope  of  the  strain-hardening  after  reverse  yielding,  and  is 
the  additional  plastic  strain  during  unloading. 

ELASTIC-PLASTIC  LOADING 

Consider  a  thick-walled  cylinder,  inner  radius  a,  and  external  radiuB  b, 
which  is  subjected  to  inner  pressure  p.  The  material  is  assumed  to  be  elastic 
Ideally  plastic,  obeying  the  Treses' s  yield  criterion  and  the  associated  flow 
theory.  The  elastic-plastic  solution  during  loading  has  been  found  by  Koiter 
(ref  7).  The  expressions  for  the  stresses  and  strains  are: 


<V0o  1 

2 

P 

(2a) 

--(51 

+  --)  -  log  -  ,  in  (a  <  r  <  P) 

Oq/Oq  2 

b2  r 

(3a) 

°r/°o  1 

2  2 

P  P 

(2b) 

m  — 

(--  5  --)  ,  in  (p  <  r  <  b) 

b2  r2 

OQ/Oq  2 

(3b) 

Vz/Oo 

-  v(or+O0)/oo  +  E Cz/a0 

(*> 

3 


(I  -2v)(14-v) 


2 


(5) 


E  u 


°o  r 


~  t  (1-v2)  -- 


and 


(E/oa)  tz  -  (,.-2v)(p/o0)/(b2/n2-l)  (6) 
where  p  »  0  (open-end),  1  (closed-end),  and  p  Is  the  elastic-plastic  boundary 
relating  to  the  Internal  pressure  p  by 


p/0o  -  -  ( l~e‘ /b  )  +  lng(n/a) 


The  equivalent  plastic  strain  can  ho  calculated  by 


(?) 


and 


(t'./u0H  -  <,  a  pi(  i>'Vr2-l)  ,  In  (a  <•  r  <  p) 

*1  =  ( 2  /  »T)  ( 1-  v-) 


(8) 


REVERSE  YIELDING 

If  the  pressure  p  given  hy  Eq .  (7)  Is  subsequently  completely  removed 
with  no  reverse  yielding,  the  unloading  Is  entirely  elastic  and  the  solution 
Is  given  by 

(9) 
(10) 


r 

°e* 


P  b 

- - !  i  - 


b2/a  -l 


oz’  -  v(or'+oo’)  +  E  uz’ 

E  tz'  =  -( p-2v)p/(b2/n2-l) 

E  u/r  =  -[(1-v-uv)  t-  (H-v)b2/r2]p/(b‘>/a‘!-l) 


(11) 

(12) 

(13) 


Let  a  double  prime  denote  a  component  In  the  residual  state,  i.e.,  ujj" 
a o  +  Oq'.  Assuming  a  reduced  compressive  yield  strength  as  a  result  of  the 
Bauschlnger  effect  and  using  Tresca’s  criterion  subject  to  or"  >  jz"  >  oq", 
the  reverse  yielding  will  not  occur  if 


r 


U 


°r”  "  °0"  -  f  °o  (U) 

Substituting  the  loading  and  unloading  solutions  Into  Eq.  (14),  we  can 
determine  the  minimum  pressure  (pm)  for  reverse  yielding  to  occur.  The 
equation  for  pm  le  given  by 

Pffl/°o  "  \  (l+f)(l-*2/b2)  (15) 

Equating  Eq.  (7)  to  Eq .  (15),  we  can  calculate  p^  and  determine  the  maximum 
amount  of  overstrain  for  reverse  yielding  not  to  occur. 

ELASTIC-PLASTIC  UNLOADING 

Now  suppose  that  the  loading  has  been  such  that  the  Internal  pressure  Is 
larger  than  pra  given  by  Eq.  (15).  On  unloading,  yielding  will  occur  for  a  < 
r  <  p'  with  p"  <  p.  Taking  into  account  the  Bauschinger  effect  (f)  and  the 
strain-hardening  during  unloading  (m'),  we  have 

or"  -  o0"  -  f  0o  +  m'Ec,P/ (l-m' )  (16) 

assuming  that  or"  >  az "  >  oq"  in  a  <  r  *  p' . 

The  material  Is  assumed  to  be  elastic-plastic,  obeying  the  Tresca's  yield 
criterion,  the  associated  flow  theory,  and  a  linear  strain-hardening  rule 
during  unloading.  Following  the  procedure  In  Bland's  work  (ref  8),  a  closed- 
form  solution  for  elastic-plastic  unloading  can  be  obtained.  The  stresses  In 
the  reverse  yielding  zone  (a  <  r  <  p' )  are  given  by 

0r'/ao  -  p/ac  -  l-  &2'  (l+f)(p'/a)2(l-a2/r2)  -  (l-32')U+f)  log(r/a)  (17) 

°e'/o0  -  V/°o  -  <1+0  "  m'cVd-m')  (18) 

d  -  6i'(l+f)(p'2/r2-l)  (19) 
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where 


^  ( 1-m' ) 

Hi'  "  (l-m’)/t®’  +  “”~2^  •  B2'  “  m’  /(1-ra’  ) 

The  stresses  in  the  elastic  tone  (p’  <  r  <  b)  are 

°r  /0°  -  -  ( 1+f ) I ±  (p’/r)2  -  (p’/b)2] 

°o’/o0  2 

The  other  expressions  for  the  entire  tube  (a  *>  r  v  b)  are 

oz'/o0  -  v(0r'+06'  )/o0  +  E  c2'/0o 
E  ez’/o0  -  -( p-2v)(p/o0)/(h2/a2-l) 


(20) 

(21) 

(22) 

(23) 

(24) 


(E/o0)  u'/r  -  (I-2u)<l+v)(or'/o0)  »  (I-v2)(l+f)(p,/r)2  -  v£  ez’/o0  (25) 
The  residual  stresses  and  the  residual  displacement  are  found  by  addition 


or  +  or’ 


09"  “  00  +  Ofl’ 
and  u”  “  u  +  u’ 


oz  +  02- 


(26) 


NUMERICAL  RESULTS  AND  DISCUSSION 

The  numerical  results  for  two  closed-end  thick-walled  cylinders  with  b/a 
■  2  and  3  were  obtained.  Various  values  of  f  and  m’  were  used  for  the  purpose 
of  showlnR  the  Bauschinger  and  hardening  effects  on  residual  stresses. 

Although  the  numerical  results  for  all  the  stresses,  strains,  and  displacement 
were  obtained,  only  those  for  the  residual  hoop  stresses  will  be  shown  here. 
The  material  constants  used  in  all  cases  were  v  -  0.3,  E/oQ  «  200.  The  slope 

of  unloading  after  reverse  yielding  was  estimated  to  be  0.3E  for  a  high 

strength  steel  (ref  14). 

The  residual  stress  distributions  in  an  autof rettaged  thick  tube  of  wall 
ratio  two  are  shown  in  Figure  3  for  p/a  ■  1.4  and  1.8.  For  a  40  percent 

autof rettaged  tube,  we  have  shown  the  results  for  three  cases  (a)  f  »  1 ,  m'  “ 
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0;  (b)  f  ■  0.50,  m'  ■  0;  (c)  f  ■  0.50,  m'  •  0.3.  The  first  case  represents  no 
Bauachinger  effect  with  no  reverse  yielding.  The  second  case  shows  the 
Bauachinger  effect  only  (ref  11).  The  third  case  shows  the  influence  of  the 
combined  Bauachinger  and  hardening  effects.  For  an  80  percent  autof rettaged 
tube,  we  have  shown  the  residual  stress  distributions  for  three  cases  (a)  f  * 
1,  m'  •  0;  (b)  f  ■  0.42,  m'  -  0;  (c)  f  ■  0.42,  m*  *  0.3.  Now  the  Influence  of 
the  Bauschinger  and  hardening  effects  is  more  significant.  Comparing  the 
residual  hoop  stress  at  the  bore  given  by  cases  (a)  and  (b)  with  case  (c) ,  the 
results  indicate  chat  negLectlng  both  effects  will  overestimate  by  46  percent, 
while  including  the  Bauschinger  effect  only  will  underestimate  by  25  percent. 

The  residual  stress  distributions  in  an  autof rettaged  thick  tube  of  wall 
ratio  three  are  shown  in  Figure  4  for  p/a  "1.4  and  2.2.  The  values  of  f  used 
for  the  20  and  60  percent  autofrettage  are  0.54  and  0.40,  respectively.  If 
the  hardening  effect  (m')  is  neglected,  we  would  have  a  smaller  compressive 
stress  at  the  bore.  If  we  neglect  both  the  Bauschinger  and  hardening  effects, 
l.e.,  f  “  1  and  m'  -  0,  we  would  have  a  larger  residual  compressive  stress  at 
the  bore.  At  60  percentage  autofrettage,  reverse  yielding  still  occurs  but  In 
a  smaller  portion  of  the  tube  around  the  bore. 

In  order  to  further  discuss  the  Bauschinger  and  hardening  effects  on  the 
residual  stress  distributions,  we  have  used  other  values  for  f  and  m'  In  a 
thick  tube  with  wall  ratio  three  and  100  percent  autofrettage.  Figure  5  shows 
the  Bauschinger  effect  (f  -  0.36,  0.68,  i.00)  only  with  no  hardening  (m'  »  0) . 
Figure  6  shows  the  effect  of  hardening  (m'  -  0,  0.1,  0.2,  0.3)  with  f  *  0.36. 
Figure  7  shows  the  Bauschinger  effect  (f  ■  0.36,  0.68,  1.00)  with  hardening 
(m'  ■  0.3).  These  results  indicate  that  the  influence  of  the  combined 
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Bauach lager  and  hardening  effects  on  the  residual  stress  distribution  is 
significant. 

CONCLUSIONS 

A  new  theoretical  model  for  a  high  strength  steel  has  been  proposed.  The 
real  Bauschinger  effect  factor  can  be  used  to  determine  the  range  of  elastic 
unloading.  The  small  strain-hardening  during  loading  Is  neglected,  but  the 
large  strain-hardening  after  reverse  yielding  is  taken  Into  account. 

A  closed-form  solution  for  calculating  the  residual  stresses  and  strains 
with  reverse  yielding  has  been  obtained.  The  numerical  results  for  the 
residual  stress  distributions  In  two  autof rettaged  thlck-walled  cylinders  have 
been  given.  The  new  results  Indicate  that  the  Influence  of  the  combined 
Bauschinger  and  hardening  effects  on  the  residual  stress  distribution  1? 
significant. 

Comparing  the  residual  hoop  stress  at  the  bore  for  an  80  percent 
autof rettaged  tube  with  wall  ratio  two,  Koiter’s  model  neglecting  both  effects 
(ref  7)  will  overestimate  by  46  percent,  while  Parker’s  model  including  the 
Bauschinger  effect  onlv  (ref  11)  will  underestimate  by  25  percent  when 
compared  with  the  present  model  taking  Into  account  both  the  Bauschlnger  and 
hardening  effects. 
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Figure  3.  Residual  stress  distribution  in  an  autofrettaged  tube 
Ib/a  =  2,  p/a  =  1.4  and  1.8). 
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Figure  6 
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Hardening  effect  on  residual  stress  distribution 
(b  =  p  =  3a,  f  =  0.36) . 
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Figure  7 
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Combined  Bauschinger  and  hardening  effects  on  residual 
stress  distribution  (b  =  p  =  3a,  m'  =  0.31. 
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